Warm dark matter at small scales: 
peculiar velocities and phase space density. 

Daniel Boyanovsky 1 'B 

1 Department of Physics and Astronomy, 
University of Pittsburgh, Pittsburgh, PA 15260 
(Dated: May 4, 2011) 

Abstract 

We study the scale and redshift dependence of the power spectra for density perturbations and 
peculiar velocities, and the evolution of a coarse grained phase space density for (WDM) particles 
that decoupled during the radiation dominated stage. The (WDM) corrections are obtained in a 
perturbative expansion valid in the range of redshifts at which N-body simulations set up initial 
conditions, and for a wide range of scales. The redshift dependence is determined by the kurtosis 
/?2 of the distribution function at decoupling. At large redshift there is an enhancement of peculiar 
velocities for fa > 1 that contributes to free streaming and leads to further suppression of the 
matter power spectrum and an enhancement of the peculiar velocity autocorrelation function at 
scales smaller than the free streaming scale. Statistical fluctuations of peculiar velocities are also 
suppressed on these scales by the same effect. In the linearized approximation, the coarse grained 
phase space density features redshift dependent (WDM) corrections from gravitational perturba- 
tions determined by the power spectrum of density perturbations and /?2- For /?2 > 25/21 it grows 
logarithmically with the scale factor as a consequence of the suppression of statistical fluctuations. 
Two specific models for WDM are studied in detail. The (WDM) corrections relax the bounds on 
the mass of the (WDM) particle candidate. 
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I. INTRODUCTION 



The current paradigm of structure formation, the ACDM standard cosmological model, 
describes large scale structure remarkably well. However, observational evidence has been 
accumulating suggesting that the cold dark matter (CDM) scenario of galaxy formation may 
have problems at small, galactic, scales. 

Large scale simulations seemingly yield an over-prediction of satellite galaxies [l[ by almost 
an order of magnitude jlH5|- Simulations within the ACDM paradigm also yield a density 
profile in virialized (DM) halos that increases monotonically towards the center [l], 049J and 
features a cusp, such as the Navarro-Frenk- White (NFW) profile [6] or more general central 
density profiles p(r) ~ r _/3 with 1 < /3 < 1.5p, @, @|. These density profiles accurately 
describe clusters of galaxies but there is an accumulating body of observational evidence jlOl 



17| | suggesting that the central regions of dark matter (DM)-dominated dwarf spheroidal 



satellite (dSphs) galaxies feature smooth cores instead of cusps as predicted by (CDM). Some 
observations suggest (l9| that the mass distribution of spiral disk galaxies can be best fit by a 
cored Burkert-type profile (l9|. This difference is known as the core-vs-cusp problem (l6l. 17 . 
The case for core-dominated halos has been recently bolstered by the analysis of rotation 
curves from the THINGS survey [F" 



Warm dark matter (WDM) particles were invoked |2CH22j as possible solutions to the 
discrepancies both in the over abundance of satellite galaxies and as a mechanism to smooth 
out the cusped density profiles predicted by (CDM) simulations into the cored profiles that 
fit the observations in (dShps). (WDM) particles feature a range of velocity dispersion 
in between the (CDM) and hot dark matter leading to free streaming scales that smooth 
out small scale features and could be consistent with core radii of the (dSphs). If the free 
streaming scale of these particles is smaller than the scale of galaxy clusters, their large scale 
structure properties are indistinguishable from (CDM) but may affect the small scale power 
spectrum 23[ providing an explanation of the smoother inner profiles of (dSphs) and fewer 
satellites. 

Furthermore recent numerical results hint to more evidence of possible small scale dis- 
crepancies with the ACDM scenario: another over-abundance problem, the "emptiness of 



voids" [24j and the spectrum of "mini- voids" [25j , both of which may be explained by a 



WDM candidate. Constraints from the luminosity function of Milky Way satellites [26] sug- 
gest a lower limit for the mass of a WDM particle of a few keV, a result consistent with 
galaxy power spectrum] 30] 



Lyman-a |27H29 



and lensing observations [3JJ. More recently, 
results from the Millenium-II simulation [32| suggest that the ACDM scenario overpredicts 
the abundance of massive > 10 10 M© haloes, which is corrected with a WDM candidate of 
m ~ 1 keV. A model independent analysis suggests that dark matter particles with a mass 
in the keV range is a suitable (WDM) candidate (H, 34]. Recent counterarguments 35, 36 
seem to suggest that (WDM) cannot explain cores in (LSB) galaxies, thus the controversy 
continues. 

In absence of conclusive evidence in favor of or against cusps or cores, and in view of 
the ongoing controversy and the body of emerging evidence in favor of (WDM), a deeper 
understanding of the small scale clustering properties of (WDM) candidates is warranted. 

Motivation and goals: 

Redshift dependence of the power spectrum and peculiar velocities: recent N-body simu 

set up initial conditions at z 



lations of WDM £51 26 



40126] or z = 50 251 with a rescaled 



version of the CDM power spectrum from a fit provided in ref. 2^] that inputs a cutoff from 



2 



free streaming, however, these simulations neglected the velocity dispersion of the WDM 
particles in the initial conditions. We seek to understand both the redshift dependence of 
the matter and peculiar velocity power spectrum in this range of redshifts for a wide range 
of scales. 

Phase space density: in a seminal article Tremaine and Gunn provided bounds on 
the mass of the DM particle from phase space density considerations: whereas in absence 
of self-gravity the fine-grained phase space density (or distribution function) is conserved 
after the DM species decouples from the plasma, phase mixing theorems [38( assert that a 
coarse grained phase space density always diminishes as a result of phase mixing (violent 



relaxation) |38|, |39 



Therefore the microscopic phase space density provides an upper bound 
from which constraints on the mass can be extracted. These arguments were generalized 
in refs. [lol. 33 . 40-43] to a coarse grained phase s pace d ensity obtained from moments of 

In ref. 1331, 



4lH43| this coarse grained phase space 



the microscopic distribution function, 
density was combined with photometric observations of (dShps) to constrain the mass and 
the number of relativistic degrees of freedom at decoupling. 

Although the microscopic phase space density, namely the distribution function, obeys 
the collisionless Boltzmann equation, the evolution of the coarse grained phase space density 
is not directly obtained from this equation (see discussion in ref.[39|). Although the proxy 



phase space density introduced in refs. (lOj, |40H42| is conserved after decoupling, its evolution 



does not include self-gravity. Therefore there remains the unexplored question of precisely 
what happens to the microscopic phase space density or its proxy introduced in refs. flU 40- 
42] when gravitational perturbations are included in the Boltzmann equation. One aspect 
is clear: the perturbations of the distribution function (microscopic phase space density) 
feature two moments that grow under gravitational perturbations: the first moment (density 
perturbations) and the second moment (velocity perturbations) which are actually related 
via the continuity equation on sub-horizon scales. In this article we study the evolution of 
the coarse-grained phase space density introduced in refs. (lol. 4oT42] as a function of redshift 
and scale for (WDM) particles in order to assess how the original arguments are modified 
by gravitational perturbations, again in the regime of redshifts at which N-body simulations 
set up initial conditions. 

Results: 

Armed with the results recently obtained in ref.jHI we obtain a perturbative expansion 
of the redshift corrections to the matter, peculiar velocity power spectra and evolution of 
a coarse-grained phase space density. This expansion is valid in the regime z <C z eq for a 
wide range of scales and is a distinct feature of (WDM) particles. These corrections depend 
on the kurtosis 02 of the unperturbed distribution function. Peculiar velocities contribute 
to the velocity dispersion and free streaming and lead to a suppression of the matter power 
spectrum for /3 2 > 1 at scales smaller than the free streaming scale at redshifts z ~ 30 — 50. 
The peculiar velocity power spectrum is enhanced at these scales and reshift leading to an 
increase of the peculiar velocity autocorrelation function and a suppression of statistical 
fluctuations. For (WDM) perturbations in the linearized approximation, it is found that 



the coarse grained phase space density introduced in refs. [10l. grows logarithmically 

with the scale factor for /3 2 > 25/21. Two specific models of (WDM) particles motivated 
by particle physics are studied in detail. Implications on the bounds for the mass of the 
(WDM) particle are discussed. 
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II. PRELIMINARIES 



We begin by establishing some notation and conventions that are used in the analysis. 
Since we focus on the region of redshift z > 1 we can safely neglect the dark energy 
component and we consider a radiation and matter dominated cosmology with 



H 2 = ^ 



Hi 



Tj2r) 

-"0 m r 

-7— [a + a, 



cq\ 



(HI) 



where the dot stands for derivative with respect to conformal time (rj), the scale factor is 
normalized to a = 1 today, and 



*eq 



Introducing 



it follows that 



1 



3229 



< ! eq 



da 
dr] 

At matter-radiation equality we define 



>'eq 



[1 



al 2 



(II.2) 
(II.3) 
(II.4) 



keq — H e q CL e g 



V2 



7T2(~) 

"0 m 



>>eq 



9.8 x 10- 
Mpc 



(II.5) 



corresponding to the comoving wavevector that enters the Hubble radius at matter-radiation 
equality, where we have used VL m h 2 = 0.134 jioT]. 

We study the evolution of perturbations in the conformal Newtonian gauge 



i.i 



goo = -a 2 {v) 1 + 2ip{x,v) 
9ij = a 2 (r]) 1 - 2cj)(x,r)) 5, 

The perturbed distribution function is given by 

f(p, X, 7]) = f (p) + Fi(p, X, T]) 



(II.6) 
(II.7) 

(II.8) 



where fo{p) is the unperturbed distribution function, which after decoupling obeys the colli- 
sionless Boltzmann equation in absence of perturbations and p, x are comoving momentum 
and coordinates respectively. As discussed in ref. (41-43] the unperturbed distribution func- 
tion is of the form 

/o(p) = /o(3/;xi,X2,--0 (n.9) 



where 



P 



(11.10) 
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where p is the comoving momentum and Tq^ is the decoupling temperature today, 



T , d ={-Y Tcmb, (11.11) 



9d 



with ga being the effective number of relativistic degrees of freedom at decoupling, Tcmb = 
2.35 x 10~ 4 eV is the temperature of the (CMB) today, and Xi are dimensionless couplings 
or ratios of mass scales. 

We neglect stress anisotropies, in which case = ip and introduce 

~ F{ pr Kn) = ^M . j ip) = m. (ma) 

n n 

where 

n ° = I Wf foiph (IL13) 

is the background density of (DM) today. Therefore 

S(k,v) = I (01^)- ( IL14 ) 

becomes Sp m /p m after the DM particle becomes non-relativistic. 

Introducing spatial Fourier transforms in terms of comoving momenta k (we keep the 
same notation for the spatial Fourier transform of perturbations), and neglecting stress 
anisotropies the linearized Boltzmann equation for perturbations is given by [46l45^| 



F(k,p;rj) + i 



k fip 
e(p,r}) 



F(k,p;T)) + 



df(p) \ 
dp J 



p <fi(k,r)) — ik pe(p,r]) <f)(k,T)) =0 (11.15) 



where dots stand for derivatives with respect to conformal time, p. = k • p, and e(p, rj) = 
\Jp 2 + m 2 a 2 (ri) is the conformal energy of the particle of mass m. The gravitational poten- 
tial is determined by Einstein's equation (461. 471 . 

As discussed in ref. [3] for a (WDM) particle with a mass in the ~ keV range, there are 
three stages of evolution: I) radiation domination and the DM particle is relativistic, II) 
radiation domination and the DM particle is non-relativistic and III) the matter dominated 
stage, during which cold and warm DM particles are non-relativistic. 

During stages (I) and (II) the gravitational potential is completely determined by the 
radiation component and the Boltzmann equation for the distribution function of the WDM 
particle is solved by integrating eqn. flII.15j) with <ft being determined by the radiation 
component. During stage (III) the gravitational potential is determined by the matter 
component and the Boltzmann equation becomes a self-consistent Vlasov-type equation. 

Since the Boltzmann equation is first order in time, the solution during stages (I) and 
(II) becomes the initial condition for the evolution during stage (III). 

In this article we focus on the evolution of peculiar velocities and phase space density 
during the matter dominated stage 10 < z < z eq , corresponding to stage (III) during which 
dark energy can be neglected. Typical N-body simulations setup initial conditions which 
input the matter power spectrum from linear perturbation theory at z ~ 30 — 50. 
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In this stage the WDM is non-relativistic, hence p/c(p, rj) = p/ma(j]), and the Bolzmann 
equation simplifies by introducing the variable 



drf __ 2V2u(r]) 



a(rf) 



where the dimensionless variable 



u (jl) = 2 ln 



y/1 + a{r}) - 1 
Jl + a( V ) + l 



u N r < u{rj) < 



(11.16) 



(11.17) 



is normalized so that w(oo) = and introduced |44l 

u NR = ln ; a NR = (V 2 (t eq ))' 



(11.18) 



where aNR corresponds to the time when the particle becomes non-relativistic, and (VHt Pq )) 
is the velocity dispersion of the DM particle at matter-radiation equality given by|44j 



(V 2 (W)^7.59 xl(T 4 V^(— ) f 1 )" 

v \ m J \g d J 



(11.19) 



In this expression g d is the number of relativistic degrees of freedom at decoupling and we 
introduced the moments 

j °° y 2+n fo(y) dy 



j^y 2 h{y)dy 

The function u[z] as function of redshift is displayed in fig. (CQ) and 

1 



a{u) 



sinh 2 [u] 



(11.20) 



(11.21) 



The solution of the Boltzmann equation during stage (III) is given in ref.([4 



F(k,p;s) 



+ im 



ds 1 a 2 (s')(j)(k,s')[k-V p f(p))e 



-i^(s-s') 



+ e" 



F{k,p; r] NR ) + (j>(k, r] NR ) ( p — 



(11.22) 
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The term F(k,p;rjNji) i n the bracket in eqn. fill. 221) is the solution of the Boltzmann 
equation at the beginning of stage (III) (end of stage (II)) its form is given in detail in 
ref.jiij but is not necessary in the discussion that follows. 

After radiation-matter equality when the WDM particle is non-relativistic and (DM) per- 
turbations dominate the gravitational potential and for k 3> k eq , the gravitational potential 
is determined by Poisson's equation 47 



#M) = ~Jm&«) (H.23) 
4 Ar a 

For s > s eq , the integral in s' in (III. 22[) is split from snr up to s eq and from s eq up to s. In 
the first integral the gravitational potential is determined by perturbations in the radiation 
fluid and in the second integral the gravitational potential is replaced by Poisson's equation 
fill. 231) . leading to the result (valid for s > s eg )jll 

+ F[k,p;s}. (11.24) 

where 

Q = —(s-s') ; /i = k-p, (11.25) 

m 

and s] is given by the second line in (111.22j) plus the contribution from the integral 

between s^r and s eq (for details see ref. |44j). 

We are interested in the corrections to the power spectra in the regime of redshift 1 < 
z <C z eq corresponding to a > 1. 

In the asymptotic limit a 1 when density perturbations grow as in an Einstein-DeSitter 
cosmology, 5 oc a, therefore in this limit and for k ^> k eq we can neglect the first term in 
( 111.24p . Since in this limit S(k,s') oc a(s') oc (1/V) 2 the integral in eqn. (1II.24|) is oc 1/s 3 
and dominates all other terms in eqn. (III. 24[) since the last term remains finite in the limit 
s ->• 00. 

Therefore in the asymptotic limit a ^> 1 and for small scales k 3> k eq the leading contri- 
bution to the perturbation in the distribution function is given by 

F(k, p; S ) ~ -i 3mk ^ f S ds' a(s') 5(k, s') ) e~^ Q (11.26) 

j s e q y 

With this form for the distribution function we can obtain any expectation value once 5(k, s) 
is determined from the solution of the Boltzmann equation. 

In ref. 44j it is shown that in terms of the variable u defined by eqns. ( 1II.16IIII.17j) 5 obeys 
the fluid-like integro-differential equation 



d 2 r u ~ 

-—5(k,u)—6a(u)5(k,u) + K 2 5(k,u) — 6a / a(u')U\a(u — u')]5(k,u') du' 
du Ju NR 

J[k,u] (11.27) 



where the inhomogeneity J[k, u] is given explicitly in ref. j44j, and 

i poo roo 

II [a(u-u')] =-y y foivW-y 2 ) sin [y a (u-u')]dy ; N = J y 2 fo(y)dy. (11.21 
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In the above expressions we introduced 



m = 2V2 * 2.15 x Itr 3 (JL) (*-)* (™) * 0.22 k (—Y (™) x (Mpc) 



VTl k e qCL e q 



and 



K 



y 2 a . 



In terms of the free streaming wavevector[44 



> eq 



2 (v 2 (t eq ))^ 



11.17 / m 



*=> VkeV 



9d\3 



(Mpc)" 1 , 



it follows that 



V6k V6X 



fs 



(11.29) 
(11.30) 

(11.31) 
(11.32) 



kf s X 

where Xf s is the free streaming length and A the wavelength of the perturbation. The (CDM) 
limit corresponds to A/ s — > 0, namely k — > 0, therefore all the WDM corrections are in terms 
of AC. 



In the (CDM) limit eqn. flll.27p reduces to Meszaro's equation 53M55| for (CDM) pertur- 
bations in a radiation and matter dominated cosmology [44j . 
The power spectrum of density perturbations is given by 



P s (k) = Ak ns T 2 (k) 



(11.33) 



where where n s = 0.963 [45j is the index of primordial scalar perturbations, A is the am- 
plitude and T(k) is the transfer function. It is convenient to normalize the WDM power 
spectrum and transfer function to CDM, namely 



Pwdm(k) — Pcdm(k) T (Ji) ] P\k) 



Pyjdmik': ^) 
Tcdm{k) 



where 



Pcdm(k) 



Ak na T 2 dm (k) 



(11.34) 



(11.35) 



is the (CDM) power spectrum and the dependence on WDM is encoded in the k dependence 
of T wdm (k; k) so that T wdm (k; k — 0) = T cdm {k). The dependence on k describes the velocity 
dispersion and non-vanishing free streaming length of the WDM particle. 

In ref . [44j it is shown that eqn. (III. 271) can be solved in a systematic Fredholm expansion, 
from which the transfer function of density perturbations at z = is extracted. The leading 
order term is a Born-type approximation which provides a remarkably accurate approxima- 
tion to the transfer function and reproduces numerical results available in the literature in 
several cases (for discussion and comparison see[44j| ) . The definition of the power spectrum 
and transfer function above are at z — 0. We seek to study the redshift dependence for 
z < 30 — 40 at which N-body simulations set up initial conditions. 

Asymptotically during the matter dominated era as a — > 00 (u — > 0) it is found 44 
that 5(k, u) — > a(u) 5(k, 0) + ■ ■ • where the dots stand for subleading terms. The leading 
and subleading asymptotic behavior in the u — > (a — > 00) limit can be obtained from 
eqn. fill. 271) . In this limit the inhomogeneity J[k, 0] is a finite constant (see expressions in 



S 



ref.[44|), the integral term receives the largest contribution for u' ~ u ~ and in this region 
we find _ 

-6«n [a(u - u')] ~ k a (u - u'f ( 1 - (3 2 ) + ■ ■ ■ . (11.36) 



where 



(11.37) 



is the kurtosis of the distribution function of the decoupled particle, with the moments 
defined by eqn. (111.20j) . and the dots stand for terms that yield subleading corrections (see 
below). 
Since 



1 



1 1 



a\u) 



- 77 + 0(u 2 



sinh 2 [u] u 2 3 
we propose the asymptotic expansion 

5(k, u) = 6 -^p- + 5i(k) + 8 2 {k) u 2 H-u\ + 



Introducing this expansion in eqn. ( III. 271) we find 

5{k,0) 



Si(k) 



[k 2 + 2] ; 5 2 (k) = -5(k,0)^-{l-P 2 ) 



(11.38) 



(11.39) 



(11.40) 



6 j -v-p-/ 4 

where 5(fc, 0) is obtained from the asymptotic solution of the full eqn. ( III. 271) . Therefore 



S{k,u) = 5(k,0)D[k,u] 
where the wavevector dependent growth factor is found to be 

D[k,u] = D cdm [u] D[k,u] 
with D c dm[u] being the CDM growth factor (for k = 0) 



D c dm n 
U 



i 11 

l + T + - 



and 



D[fc,u] 



(km) 2 (km) 4 . , r 1X , 
l + l_^ + l_^(-ln[-M]) (l-/3 2 ) +■ 



(11.41) 



(11.42) 



(11.43) 



(11.44) 



contains the WDM corrections as is manifest in the k dependence. 
For m — y we find 

11.2 



D cd m[u] 



— a + 
2 3 3 



which is recognized as the growing solution of Meszaros equation for CDM 53H55|. Further- 



more from (III. 161) we recognize that 



— K U 



kl fs[ ,Vo,V\ 



m 



(11.45) 
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r a / it) ) "i 

where lf s [ m iVoiV] is the comoving free streaming distance that a particle with (comov- 
ing) velocity a/ (p 2 }/m travels between conformal time r] and today 770 3> 1. We see that 
up to logarithms, the expansion in powers of k u is valid at late times for wavelengths much 
larger than the free streaming distance that the particle would travel between that time and 
today. 

The identification (1II.45P leads to a simple physical interpretation of the first term in 
D[k,u\: free streaming of collisionless particles suppresses the gravitational collapse of den- 
sity perturbations, the longer the time scale, the farther the free streaming particles can 
travel away from the collapsing region erasing the perturbations. Therefore the first term 
reflects that at earlier times (larger values of u) density perturbations are larger. The sec- 
ond term, however, has a more interesting interpretation. As it will be discussed below, 
it represents the peculiar velocity contribution to free streaming induced by gravitational 
self- interaction (see discussion on peculiar velocity below). When /3 2 > 1 the peculiar ve- 
locity contribution increases the free streaming velocity leading to a suppression of power, 
which counterbalances the enhancement by the first term. Which term dominates depends 
on the scale k, the free streaming wavevector, a characteristic of the WDM particle, and the 
redshift. This will be analyzed in two specific models below. 

We emphasize that the expansion in ( 1II.39j) is valid at long time, in particular for km < 1. 
At higher orders in the expansion, the terms that feature the ln(— u) only appear linearly 
in the logarithm but multiplied by higher powers of ku, therefore for \ku\ < 1 the third 
term in D is the leading logarithmic contribution, with higher contributions being of the 
form (nu) n ln(— u) ; n = 6, 8 • • • . This is an important observation: in particular within the 
regime of validity of the perturbative expansion \nu\ < 1 it is still possible that \ku ln(— u) \ ~ 
1 and the second term in (1II.44I) can balance the first term within the region of validity of 
the approximation. 

An estimate of the range of validity is obtained from 

0.098 < -u[z] < 0.125 for 30 < z < 50 , (11.46) 

for example in the region of redshifts where initial conditions for N-body simulations are 
set up, the WDM corrections to the growth factor are of O(10 — 15%) for k > (1 — 2) kf s 
which for a species with m ~ keV decoupled with g d ~ 30 — 100 with y 2 ~ 10 corresponds 
to k > 10 - 30 (MpcT 1 . 

There is a caveat in this analysis of the reliability of the expansion, since it applies only in 
the linear regime where the linearized Boltzmann equation describes the transfer function. 
It is conceivable that non-linear effects restrict further the regime of validity, but of course 
this cannot be assessed in the linear theory which is the focus of this discussion. 

Using Poisson's equation (III.23D . the asymptotic behavior 5(k,u) — > 6(k,0) a{u) and the 
definition of the transfer function [471] T(k) 

0(M >> 1) = y^)T(^) (11.47) 

where 4>i(k) is the primordial value of gravitational perturbations seeded by inflation. It 
then follows that 

5(k,0) = -Uk)^T(k). (11.48) 
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We emphasize that there are two different averages: i) the statistical average of a quantity 
O with the perturbed distribution function fo + F\ to which we refer as (O), ii) average over 
the initial gravitational potential 0j which is a stochastic Gaussian field (we neglect possible 
non-Gaussianity) whose power spectrum is determined during the inflationary era 

<j>i(k)<t>i{-k') = (2tt) 3 6® (k - k') P^k) (11.49) 

where the A B refers to averages with the primordial Gaussian distribution function for the 
gravitational potential 1 . Therefore full expectation values correspond to averages both with 
the perturbed distribution function and the Gaussian distribution function for the primordial 
gravitational potential, these are given by (O), with the power spectrum of matter density 
fluctuations 

8(k,0)S(-k',0) = (2ir) 3 5®(k- k')P s (k), (11.50) 

where Ps(k) is given by eqn. fill. 331) . 

Including the wavevector dependent growth factor D[k,u] ( 1II.44I) but keeping only the 
(WDM) (k / 0) corrections with redshift, the effective (WDM) power spectrum at z <C z eg 
is given by 

P s [k,z] = P wdm (k)D 2 [k,z] (11.51) 
where the scale and redshift dependent correction is given by (see eqn. (III.44I) ) 



D[k,z] 



K 2 


1 + z 


K 4 


1 + z 


2 

In 


1 + Z eq 


~6~ 


1 + z eq 


8 


1 + Z eq 


1 + z 



- 1) + 



(11.52) 



For /3 2 > 1 the third term is negative and competes with the second term, dominating 
the corrections for scales 



k > k 



fs 



2(1 + z, 



eq ) 



9(l + z)(/3 2 -l)ln 



(l + Zeq) 
(1+*) 



(11.53) 



for 02 — 1~0(1) and z ~ 30 — 50 one finds that the third term dominates over the second 
for k ~ (1 — 2) kf s . These are the scales beyond which the contribution from the peculiar 
velocities to free streaming leads to a suppression of the power spectrum. Coincidentally this 
is the scale at which the power spectrum displays (WDM) acoustic oscillations which arise 
from the competition between free streaming and gravitational collapse in the collisionless 
regime as described in ref. 44j . 



1 This definition should not be confused with that of the moments in eqn. (jll. 20[) which refer to averages 
with the unperturbed distribution function. The meaning of averages is unambiguously inferred from the 
context. 



11 



III. PECULIAR VELOCITY AND PHASE SPACE DENSITY: 



Statistical averages of observables with the perturbed distribution function ( III. 8 [) in the 
linearized theory (in terms of their spatial Fourier transform) are given by 



0{k;r,) = (0(p,k,r])) 



r d 3 p 

J (27T) :i 


fo(p) + Fi(p,k;r]) 


0(p, k; rj) 


O (k- V ) + AO(k- V ) 


I 


(2tt) 3 


fo(p) +F 1 (f,k; 77) 




1 + 





where 



O (k;v) 
AO{k;r]) 



d 3 p 
(2^)3 

d 3 p 
(27p 



f(p)0(p,k;r)), 
F(p, k; rj) 0(p, k; 77) . 



(HI 1) 

(III.2) 
(III.3) 



where F, f are defined in eqn. (1II.12j) . In the linearized approximation 



0{k- 77) ~ O (k; V ) + [AO(k; V ) - O (k; rj) 5{k, v) 



(III.4) 



With F(k,p;s) given by ( III. 261) and S(k,r]) given by ( lll.39jll.40l) we can now obtain any 
statistical average by expanding 0(p,k;r]) = 0{p,k, n;rj) in Legendre polynomials in /1 
and carrying out the integrals in p, \i leading to an expansion in spherical Bessel functions. 
However, here we focus on obtaining the leading asymptotic expansion of these averages for 
z z eq , namely for u 1. This is readily achieved by using the asymptotic expansion 



( 1H.39|) with the coefficients given by ( 1H.40|) . expanding 



1 i 

exp[-ifiQ] ~ 1 - ifiQ - -fi 2 Q 2 + -/i 3 Q 3 + 

2 o 

and integrating over fi and p term by term in the expansion. 



A. Peculiar velocity 

Writing the comoving peculiar velocity in terms of the longitudinal and transverse com- 
ponents 

0(k,rj) = (—) = VT + kv L , k-v T = (III.5) 
m 

where 

Vl = — a ; u = k-p (III. 6) 

m 

and p is the comoving momentum. In the linearized approximation, the expectation value 
of k vl is given by 

* Mk, v) = J tSs%> ^; v) ■ (IIL7) 
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Furthermore, F(p, k; rj) is a function of k and k ■ p leading to vt = in the linearized 
approximation. Since the gravitational potential is only a function of k, the first term on 
the right hand side of ( III. 221) does not contribute and we find 



kv L (k,T]) = i 



d 



d?p 



F(p,k;s) + ( f>(k,s){p^-) 



ds J (2tt) 3 

Using d/ds = ad/dr] equation (1111.81) becomes 



. d 
ds 



S(k, s) — 3(f)(k, s) 



(III.8) 



(III.9) 



dS d<p .k 

3— + i- v L = 

dr] arj a 

which is recognized as the continuity equation in presence of the gravitational potential (47| 2 
for the comoving longitudinal velocity. For a ^> 1 and k ^> k eq the second term in the 
continuity eqn. ( 1III.9I) can be safely neglected, leading to 



v L (k,u) = i 
As a function of redshift we find 

v L (k, z) = i 



k eq a eq dS(k,u) 
2y/2k du 

k eq a eq 5(k,0) 



y/2k (-«N) 3 

where we used the asymptotic expansion (III.39|II.40I) and introduced 



(111.10) 



(ni.li; 



K 

1 + -T 



1 + Z ]2 



l + Z, 



eq 



In 



l + z, 



eq 



l + Z 



02 



1 +■■ 



(111.12) 



In the CDM limit k — > the growth factor u 



-3 



which is recognized as the growth 



of comoving peculiar velocity in a matter dominated cosmology, in this limit T(k) — >■ 1(4 
and V w d m [k, u] — > 1. The function V w dm[k, u] encodes the corrections to the peculiar velocity 
at small scales. It is clear that as compared to the CDM case, when the kurtosis 02 > 1 the 
peculiar velocity at small scales k > 1 is larger at higher redshift. Comparing eqn. (IIII.12j) 
with the third term in eqn. (III.52j) confirms the interpretation of the suppression of the 
power spectrum at small scales and high redshift as a consequence of the peculiar velocity 
contribution to free streaming. 



B. Statistical fluctuations and correlation functions: 

In the linearized approximation (and with adiabatic perturbations only), the perturbation 
in the distribution function F(k,p; u) is linear in the primordial gravitational potential <pi{k) 
which is a Gaussian variable determined by the power spectrum of perturbations during the 
inflationary stage (here we neglect possible non-gaussianities). Therefore as discussed in the 



Note that the Newtonian potential in eqn. (|II.7[> features a minus sign with respect to the definition in 
47|. 
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previous section there are two different averages, i) a statistical average with the perturbed 
distribution function /o + F\, ii) with the initial Gaussian probability distribution of P^{k) 
in eqn. (Ill .49ft . 

Statistical fluctuations are contained in the variance of the various quantities calculated 
with the perturbed distribution F. These are linear in S(k, 0), namely linear in 0,, therefore 
they feature Gaussian fluctuations with the probability distribution function P<p(k), but with 
non-gaussian statistical variances. 

As an example of a statistical fluctuation consider 



p z 



A< = 6(i-) / du'a(u')5(k,u') 



K 



(u-u')-'-p 2 {u-u'f 
o 



du' 



(111.13) 



Using the asymptotic expansions f lII.38IIII.39pjH.40p we find up to leading logarithmic order 

5{k,0) 



a(u') 5(k, u')(u — v!) du' 



6 u 2 



1 — k u ln[— u] 



+ 



d(u')5(k,u')(u-u') 3 du = -S(k,0) \n[-u] + 



(III. 14) 
(111.15) 



leading to similar statistical fluctuations for the total velocity dispersion (p 2 /m 2 ) and the 
transverse component vt, namely (all quantities are comoving) 



A<«2> 



(p 2 )o S(k,0) 



1 - k 2 u 2 hx[-u] (l - 



A- 



,p_ 

v 9 



5 /P 2 S(k,0) 



3 m 2 



1 - k 2 u 2 \n[-u}(l 



2 (p 2 ) S(k, 0) 



3 m 2 



21 

25 

3 



+ ■ 



1-k 2 u 2 ln[-u](l - ~# 



+ 



(111.16) 



(III. 17) 



(111.18) 



Restoring units, writing (p 2 )o = y 2 T^ d and expressing these expressions in terms of redshift, 
we find the following statistical fluctuations 



A<»2> = 16.33 (i^) 2 V*( M ) 
Vsec/ V m / 



1 + z, 



(■q 



l + z 



1-K 2 



l + Z 



l + z, 



eq 



hi 



l + Z, 



eq 



l + Z 



(iii.i9) 



\ I 



Vsec/ V m / 



l + z, 



eq 



l + Z 



21 2 

1 fir 

25 



■ l + z - 

.1 + Z e(J . 



hi 



1 + z eq 

.l + Z . 



(111.20) 



25 ^ 



+ • 



A(v 2 ) 



10.89 (^) 2 (^)V*(M) 
Vsec/ \ m J 



l + z, 



L l + z J 



1 3 2 
1 K 

5 



l + z 



l + z, 



eq 



In 



l + Z J 



5 

! ~3 
(III. 2 
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These expressions also show that the WDM corrections (proportional to k 2 ) suppress the 
statistical fluctuations at small scales k > kf s where the peculiar velocity contribution to 
free streaming becomes important (we will see below that at least for the (WDM) candidates 
considered here f3 2 > 2). 

The peculiar velocity autocorrelation function is given by 



d 3 k 



Ak-x 



(2tt)< 



d 3 k' 



-ik'-x' 



Vi(k,u) v*(-k',u) 



using (111501) and (10.101) we find 

d 3 k 



P n 2 
2u 6 



p ik ' f ' h h 



\k z) 



(2tt) 3 1 J k 2 
Since there is only one vector r we write 

Zv(f, = ^(f) ^(r; z) + Vjffl e«(r; z) 



r = x — x 



where 



n Tj . 



are the projectors on directions parallel and perpendicular to r. 
We find 



k 2 n 2 
^eq^eq 



12 7T 2 (u[z]f 

k 2 n 2 

^eq^eq 

6tt 2 {u[z\Y 



dkPs(k) V wd m[k } z) jo(kr) - 2j 2 (kr 
dkP 5 (k) V wdm [k, z] j (kr) + j 2 {kr 



(111.22) 



(111.23) 

(111.24) 
(111.25) 



(111.26) 
(111.27) 



where jo,2 are spherical Bessel functions. 

Thus we see that the effective (WDM) power spectrum for peculiar velocities is 
P 5 (k)V wdm [k,z}. 

From expression 12[) it is clear that for (3 2 > 1 (WDM) perturbations enhance the 
peculiar velocity autocorrelation function for z ~ 30 — 50. This enhancement of the velocity 
correlation function is in concordance with the suppression of the power spectrum, since the 
larger velocity dispersion induced by self-gravity leads to a larger free streaming velocity 
and a further suppression of the power spectrum. 



C. Phase space density: 



In their seminal article Tremaine and Gunn [37J argued that the coarse grained phase 
space density is always smaller than or equal to the maximum of the (fine grained) micro- 
scopic phase space density, namely, the distribution function, allowing to establish bounds 
on the mass of the DM particle. 



Such argument relies on a theorem 



, |39[ that states that collisionless phase mixing 
or violent relaxation by gravitational dynamics (mergers or accretion) can only diminish 



33 



the coarse grained phase space density. A similar argument was presented in refs. [10| _ 
III! where a proxy for a coarse-grained phase space density in absence of gravitational 
perturbations was introduced. 
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However, whereas the distribution function obeys the collisionless Boltzmann equation, 
Dehnen 39) clarifies that the coarse grained phase space density does not necessarily evolve 
with the collisionless Boltzmann equation, and introduces an excess mass function which is 
argued to always diminish upon gravitational phase space mixing. 

Numerical simulations confirm the evolution of a coarse grained phase space density 
towards smaller values during violent relaxation events such as encounters, mergers and 



accretion of haloes [56l. |57||. In the simulations in ref.[56j a phase space density Q is obtained 
by averaging p, a over a determined volume, and its evolution with redshift is followed from 
z — 10 until z = diminishing by a factor ~ 40 during this interval. 

However, to the best of our knowledge, a consistent study of the evolution of the micro- 
scopic phase space density including gravitational effects even in the linearized approxima- 
tion has not yet been provided. 

In linearized theory, the corrections to the distribution function Fi, or rather the nor- 
malized perturbation F(p, k, rj) defined by eqn. ( 1II.12j) obeys the collisionless Boltzmann 
equation ( 1H.15|) . whose solution in the regime when the DM is non-relativistic is given by 
eqn. ( 1II.22|) . Thus the time evolution of the microscopic phase space density is completely 
determined. Two aspects of this solution invite further scrutiny: i) density perturbations 
grow from self-gravity effects, ii) peculiar velocities also grow, a direct consequence of the 
continuity equation (1III.9P and explicitly shown by eqn. ( IIII.lip . That both quantities grow 
upon gravitational collapse suggests an examination of the phase space evolution in the 
linearized regime. 

In principle one could perform the Fourier transform back to (comoving) spatial co- 
ordinates and obtain F{p,f;rj), however, S(k,rj) is a stochastic variable with a Gaussian 
probability distribution determined by the power spectrum of the primordial gravitational 
potential. Therefore, the linear correction to the microscopic phase space density itself 
becomes a stochastic variable as discussed above. 

Rather than pursuing the Fourier transform, which in the linearized approximation can 
be performed at any state in the calculation, we follow refs. (Tol. 33 , 40-43] and define the 
coarse grained (dimensionless) primordial phase space density 



V 



n(t) 



(111.28) 



where Pj = p/a(t) is the physical momentum. In absence of gravitational perturbations, the 
(unperturbed) distribution function of the decoupled species is frozen and n(t) = no/a 3 (t), 



therefore it is clear that D is a constant, namely a Liouville invariant, 
gravity it is given by 



In absence of self- 



Vn 



f™y 2 fo(y)dy 



2tt 2 



J™y 4 fo(y)dy 



(111.29) 



where fo(y) is the decoupled distribution function, and g the number of internal degrees of 
freedom of the WDM particle. 
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When the particle becomes non-relativistic pit) = m n(t) and (V 2 ^ = (-^), therefore, 

P Qdh 



V 



ITT 



'y 2 Y 



rrr 



(111.30) 



where Qdh = p/(V 2 ) 2 is the phase-space density introduced in refs. [Io|, 40 

In the non-relativistic regime V is related to the coarse grained phase space density Qtg 
introduced by Tremaine and Gunn 37 



Q 



TG 



P 



m 4 (2 7TC7 2 )3 



— V 

2tt 



(111.31) 



where a is the one-dimensional velocity dispersion. The observationally accessible quantity 
is the phase space density p/cr 3 , therefore, using p = ran for a decoupled particle that is 
non-relativistic today and eq.f lHI.30]) . we define the primordial phase space density 3 



Pdm 



32 m A V = 6.611 x 10 8 V 



a 



DM 



m 



keV 



4 M /kpc E 
(km/s) 3 



(111.32) 



In refs. |33l. l4lH43l| the phase mixing theorem was invoked to argue that the observed phase 
space density is smaller than the primordial value (lIII.32j) with D replaced by P given by 
eqn. (IIII.29p . lea ding to lower bound on the mass of the (WDM) particle. However, as 
emphasized in ref. 39] the phase mixing theorem 38j does not directly address the evolution 
of T>, nor has there yet been an analysis of its evolution in the linearized regime. 

The results obtained above allow us to directly calculate the corrections to T> from self- 
gravity in the linearized theory. Using the identity (1III.3j) for linearized statistical averages, 
it is given by 



V 



n 



1 + 6(k, u) 



( P 2) + A(^) 



V 



21 
20 



5(k,0)K 2 In 



l + z, 



l + z 



25 
21 



(111.33) 



were we have used eqns. ( llL39|II.40j) and ( 1III.17j) . In the (CDM) limit k — > this coarse 
grained phase space density remains constant at least up to linear order in gravitational 
perturbations. However, for (WDM) eqn. (1III.33j) clearly indicates that in the regions 
where matter density perturbations are positive the phase space density increases with the 
logarithm of the scale factor when f3 2 > 119. We will see below that this the case at least 
for two examples of (WDM) candidates supported by particle physics models. The reason 
for the increase in the coarse grained phase space density can be tracked to the suppression 
of statistical fluctuations: the leading term in S(k, u) = 5(k,0)/u 2 + ■■■ cancels against 
the leading term proportional to 5(k,0)/u 2 in the statistical fluctuation 1 7j) . these are 
the only contributions that remain in the (CDM) limit, however the (WDM) contribution 



keV 4 (km/s) 



1.2723 10 8 

kpc 
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suppresses the statistical fluctuation of the velocity dispersion leading to an increase of the 
coarse grained phase space density as a consequence of the suppression of the statistical 
fluctuations (statistical variance) of the velocity dispersion in the (WDM) case. 



IV. TWO SPECIFIC EXAMPLES. 



We now focus on two specific examples of WDM candidates: sterile neutrinos produced 
via the Dodelson Widrow (DW) mechanism j58|] for which 



fdw{y) 



ey + 1 



(IV. 1) 



where the constant £ is a function of the active-sterile mixing angle 58], and sterile neutri- 
nos produced near the electroweak scale via the decay of scalar or vector bosons (BD) for 
which El, |E3 



fbd(y) 



n=l 



-ny 



5 

712 



(IV.2) 



where A ~ 10~ 2 . 

We implement the Born approximation to the matter power spectrum presented in ref. 44 
to obtain the corrected power spectrum normalized to (CDM) [T(k)D[k, z]] 2 . As discussed 
in ref.jHJ the Born approximation yields excellent agreement with the power spectrum 
obtained in ref. 29j for (DW) sterile neutrinos. 



(DW) sterile neutrinos: 

For the distribution function (1IV.1I) we find: 



y 2 = 12.939 ; fa = 2.367 ; k 



fs 



5.44 



/ m 



9d 



keV/ V 10.75 



Mpc 



-i 



(IV.3) 



The (DW) case is displayed in fig.©: the fig. for [D[k, z]] 2 for the "standard" value 
gd = 10. 75 [HI clearly shows the crossover from an early enhancement to a later suppression 
of the power spectrum as a consequence of the contribution from peculiar velocity at small 
scales. For m = IkeV (the value used in the figure) kf s = 5.44 (Mpc) -1 , and the figure 
clearly shows that the crossover from enhancement to suppression occurs at k « 1 — 2kf s 
for 30 < z < 50 . The corrections from D[k, z] are not resolved in the log- log scale, however 
a linear-linear display of the region k > 2 kf s reveals the 10 — 15 % suppression of the power 
spectrum. This range of small scales is where the power spectrum develops the oscillatory 
behavior associated with the (WDM) acoustic oscillations discussed in ref. (44| . 

(BD) sterile neutrinos: 

Sterile neutrinos produced by the decay of scalar or vector bosons at the electroweak 
scale|43|, |59[ are colder for two reasons, i) their decoupling occurs when g d ~ 100 and they 
do not reheat when the entropy from other degrees of freedom is given off to the thermal 
plasma, ii) their distribution function (IIV.2[) is more enhanced at small momentum thereby 
yielding smaller velocity dispersion. For this species 



y 2 = 8.509 ; 2 = 2.890 ; k fs = 14.107 



keVJ V100/ 



Mpc 



(IV.4) 
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DW, m=l keV, g d =10.75 



DW, gd=10.75 , m=l keV 




10 12 

k (Mpc) 1 



DW, m=l keV, g d =10.75 



k(Mpc) 

DW, m= 1 keV, g d =10.75 




13 14 

k (Mpc) 



k (Mpc) 



FIG. 2: DW, m = 1 keV,g d = 10.75: upper left panel D[k,z] and [T(k)D(k, z)] 2 , lower panel: 
small scale region of [T(k)D(k, z)} 2 and V wdm [k, z] all for z = 30, 40, 50. 



This case is displayed in fig.©: the fig. for [/}[&;, z]] 2 for g d = 100 (corresponding to 
freeze-out at the electroweak scale) also shows the crossover from an early enhancement as a 
consequence of free streaming to a later suppression of the power spectrum as a consequence 
of the extra contribution to free streaming from peculiar velocity at small scales. For m = 
IkeV (the value used in the figure) kf s = 14.107 (Mpc) -1 , and the figure clearly shows 
that, again, the crossover from enhancement to suppression occurs at k ~ 1 — 2kf s for 
30 < z < 50 . The corrections from D[k,z] are not resolved in the log-log scale of the 
power spectrum, however a linear-linear display of the region k > 2 kf s reveals the 10 — 15 % 
suppression of the power spectrum. In this region the figure displays a hint of the (WDM) 
acoustic oscillations discussed in ref. jUj]. As discussed in ref. 44] the smaller amplitude of 
the (WDM) acoustic oscillations as compared to the (DW) case are a reflection of the fact 
that (BD) sterile neutrinos are colder as explained above. 

In both these cases, we see that there is a suppression of the power spectrum for z ~ 30—50 
in the small scale region k ~ (1 — 2) kf s and an enhancement of the peculiar velocity in the 
same region, both effects are at the 10 — 15% level and clearly correlated: the larger peculiar 
velocity adds to free streaming depressing the power spectrum. Although these effects are 
at the level of few percent, it is conceivable that they may be magnified by the inherent non- 



19 




FIG. 3: BD, m = 1 keV,^ = 100: upper left panel D[k, z] and [T(k)D(k, z)] 2 , lower panel: small 
scale region of [T(k)D(k, z)} 2 and V wdm [k, z] all for z = 30, 40, 50. 



linearities in the process of gravitational collapse, perhaps leading to important consequences 
for galaxy formation in N-body simulations. 



V. CONCLUSIONS 



Motivated by recent and forthcoming N-body simulations of galaxy formation in (WDM) 
scenarios, we set out to study the redshift corrections to the matter and peculiar velocity 
power spectra and corrections to the phase space density from gravitational perturbations 
in the region 30 < z < 50. This is the region in redshift where N-body simulations set up 
initial conditions and the dark energy component can be safely neglected. 

Drawing from results in ref . 44j , we implemented a perturbative expansion for the redshift 
and scale dependence of the distribution function, matter density perturbations and coarse 
grained phase space density valid for z/z eq <C 1 and a wide range of scales, up to leading 
logarithmic order in the scale factor. 

We find that for (WDM) the redshift dependence is determined by fa, the kurtosis of the 
unperturbed distribution function after freeze-out, with an enhancement of of the peculiar 
velocity power spectrum and autocorrelation function at larger redshift for fa > 1. This en- 
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hancement in the peculiar velocity hastens free streaming and leads to a further suppression 
of the matter power spectrum for k > (1 — 2) kf s , where kf s is the free streaming wavevector. 
For (WDM) gravitational perturbations lead to a suppression of the statistical fluctuations 
of velocities when /3 2 > 5/3. 

We also study the linear corrections to the coarse grained phase space density introduced 
m refs. HE E3-S3 

resulting from gravitational perturbations. We find that whereas these 
vanish for (CDM) resulting in a constant (coarse grained) phase space density, (WDM) per- 
turbations lead to a logarithmic growth with scale factor as a consequence of the suppression 
of statistical fluctuations if fii > 1.19. 

Two specific examples of (WDM) candidates are studied in detail: sterile neutrinos pro- 



duced non- resonantly either via the Dodelson-Widrow mechanism |58j or via the decay of 



scalar or vector bosons at the electroweak scale [43j, [59|. In these cases we find that the 
corrections to the power spectra of matter and peculiar velocities are of order 10 — 15 % for 
scales k ~ (1 — 2) kf s and redshifts z ~ 30 — 50. 

Impact on the bounds on the mass: The scale and redshift dependence of 
the power spectra are encoded in the effective matter and velocity power spectra 
P s {k)D 2 [k,z] ; P s (k)V[k,z] with D[k, z] ; V[k, z] given by eqns. ( 111.5211111. 121) respectively. 

To assess the impact of the above results on the bounds on the mass of the (WDM) 
particle consider two N-body simulations with a particle of the same mass both setting up 
initial conditions at the same z ~ 30 — 50, one with the matter and peculiar velocity power 
spectra at z ~ 1 and the other with the spectra corrected by the scale and redshift dependent 
factors obtained above. If > 1 the corrected matter power spectrum features the (WDM) 
suppression and the peculiar velocity power spectrum features the (WDM) enhancement 
for k > kf s found above. These effects at small scales are akin to the suppression of density 
fluctuations and enhancement of velocity dispersion associated with a lighter particle for 
the un-corrected power spectra. This is because a lighter particle features a smaller kf s 
and a larger velocity dispersion. Therefore these corrections allow larger masses to describe 
the same large scale structure output from the N-body simulations as compared to the 
un-corrected power spectra. Thus one aspect of the corrections is to allow larger mass 
(WDM) particles, thereby relaxing the bound on the mass, at least for those models for 
which /?2 > 1. However, this is not all there is to the corrections, because the coarse-grained 
phase space density increases, which would correspond to a colder particle with smaller 
velocity dispersion. Thus the net effect of the corrections cannot be simply characterized as 
being described by an increase or decrease of the mass of the particle and ultimately must 
be understood via a full N-body simulation. 

Although these corrections are relatively small, non-linearities arising from gravitational 
collapse may result in a substantial amplification of these effects, if this is the case, and only 
large scale N-body simulations with the corrected power spectra can assess this possibility, 
then it is conceivable (and expected) that the bounds on the mass of the (WDM) particle 
may need substantial revision. 

The results obtained here suggest a breakdown of perturbation theory either at large 
redshift and or small scales k ^> kf s , this is clearly an artifact of the expansion, the integral 
in (1II.27P which yields the logarithmic contribution is bounded and well behaved both in 
the small scale and u — > u^r limits 44|. However a systematic study of smaller scales and 
or larger redshifts would require a full numerical solution of the integro-differential equation 
( 1II.27j) . If future N-body simulations find that the corrections obtained here do modify the 
dynamics of large scale structure formation in (WDM) models substantially, such a study 
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may be worthy of consideration. 
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